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ABSTRACT 

We present a new determination of the optical luminosity function (OLF) of active 
galactic nuclei (AGN) at low redshifts (z < 0.15) based on Hubble Space Telescope 
(HST) observations of X-ray-selected AGN. The HST observations have allowed us to 
derive a true nuclear luminosity function for these AGN. The resulting OLF illustrates 
a two-power-law form similar to that derived for QSOs at higher redshifts. At bright 
magnitudes, Mb < —20, the OLF derived here exhibits good agreement with that 
derived from the Hamburg/ESO QSO survey. However, the single power law form for 
the OLF derived from the Hamburg/ESO survey is strongly ruled out by our data 
at Mb > —20. Although the estimate of the OLF is best-fit by a power law slope 
at Mb < —20.5 that is flatter than the slope of the OLF derived at z > 0.35, the 
binned estimate of the low redshift OLF is still consistent with an extrapolation of 
the z > 0.35 OLF based on pure luminosity evolution. 



1 INTRODUCTION 



The QSO optical luminosity function (OLF) and its evolu- 
tion with redshift has been studied extensively for over three 
decades (see e.g. Schmidt 1968, Marshall et al. 1983, Boyle 
et al. 1988, Hewett, Foltz & Chaffee 1993, La Franca & Cris- 
tiani 1997). This has led to a detailed picture of the QSO 
OLF over a wide range in redshift from z ~ 0.3 to z > 4. 
In contrast, the local (z < 0.15) QSO OLF is actually much 
more poorly determined, frustrating attempts to link QSO 
evolution at moderate to high redshifts with nuclear activity 
in galaxies at the present epoch. 

This is due to a number of factors associated with the 
compilation of a suitable sample of local active galactic nu- 
clei (AGN) with which to derive the local OLF. First, lo- 
cal AGN are relatively rare. Their space density is approxi- 
mately 100 times less than that of normal galaxies, and large 
area surveys are required to yield a statistically useful sam- 
ple. Secondly, many selection techniques for local AGN suffer 
from morphological biases. While surveys for stellar-like ob- 
jects are clearly biased against resolved AGN, galaxy-based 
surveys are equally biased against objects with a dominant 
nuclear component. Finally, accurate knowledge of the nu- 
clear OLF requires accurate subtraction of the light from the 
host galaxy. In the low luminosity AGN (Mb > —23) that 
constitute the vast majority of the low redshift population, 
the light from the host galaxy may dominate the nuclear 
luminosity. Even at relatively low redshifts (z ~ 0.1), see- 
ing limitations imposed by ground based observations limit 
accurate modelling of the luminosity profiles of the central 
regions of AGN host galaxies to scales typically larger than 
lh^kpc. 

A recent attempt to estimate the local AGN OLF has 



been carried out by Kohler et al. (1997), hereinafter K97. 
Using a sample of 27 candidates selected from the Ham- 
burg/ESO objective prism survey, K97 derived a local AGN 
LF that exhibited a featureless power law form over a wide 
range in absolute magnitudes —24 < Mb < —18. The form 
of the low redshift AGN LF is thus very different from 
the two-power-law luminosity function at higher redshifts 
(z > 0.5). This is a significant challenge for any theoretical 
model which seeks to connect the evolution of QSOs at high 
redshift with the local AGN population. 

The Hamburg/ESO survey covers an extensive area 
(611deg 2 ; now extended to 3700deg 2 , see Wisotzki 2000) and 
is free of morphological bias. Unfortunately the spatial res- 
olution (1-2 arcsec) of the survey is not sufficiently good to 
permit an accurate deconvolution of the galaxy and nuclear 
light even for the lowest redshift AGN (z < 0.1) in the sam- 
ple. For the 0.07 < z < 0.3 sample K97 used small-aperture, 
zero-point corrected B band CCD magnitudes which were 
subsequently corrected to reflect nuclear luminosities by sub- 
tracting a template host galaxy value of Mb = — 21 ; for the 
AGN with z < 0.07 corrections were calculated individually 
and ranged from 0.21 to 1.61 mag. 

Until recently, AGN data sets studied with HST were 
either too small or the samples on which they were based 
were too heterogeneous to construct a reliable estimate of 
the local OLF. We report here on the estimate of the local 
nuclear OLF based on HST observations of 76 AGN selected 
from a unbiased sample of X-ray selected AGN. The sample 
is part of the extensive Einstein Medium Sensitivity Sur- 
vey (EMSS, Stocke et al. 1991) which covers over 400deg 2 
and has near-complete (> 96 per cent) optical spectroscopic 
identification with no morphological bias. An earlier attempt 
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to derive the low redshift OLF using the EMSS was made 
by della Ceca et al. (1996). They used 226 broad line AGN 
with z < 0.3 to obtain a total (nuclear + host) OLF. This 
OLF was then convolved with the observed distribution of 
nuclear-to-total flux ratios for Seyfert 1 and 1.5 galaxies to 
yield a nuclear OLF. In this paper we propose to improve 
significantly on this work by using the HST observtions to 
correct explicity for the host galaxy light in each AGN. In 
section 2 we report on the measurement of the OLF from 
this sample, and in sections 3 and 4 we present and discuss 
the results obtained, comparing them to the K97 results. We 
present our conclusions in section 5. 



2 ANALYSIS 
2.1 The data 

Details of the comprehensive HST imaging survey from 
which the sample of AGN used in this analysis were drawn 
is presented in a paper by Schade et al. (2000, hereinafter 
SBL). A full discussion of the methods used to select and 
observe these AGN is presented by SBL, thus only a brief 
description will be given here. 

HST observations of 76 z < 0.15 AGN selected from 
the EMSS survey were carried out in the F814W (/) band, 
chosen to assist in the detection of the redder host galaxy 
components over the bluer nucleus. These data were com- 
plemented by deeper ground-based observations in the B 
and R bands for 69 AGN in the survey. A simultaneous 
three-component parametric model fit to the B, R and I 
images was performed for each AGN in the sample to de- 
rive magnitudes for the nuclear point source, bulge and disk 
components in each object. Despite the improved spatial res- 
olution afforded with the HST, the fitting procedure is com- 
plex and error estimation required significant modelling of 
the fitting process. For host-dominated objects uncertainties 
in Mb (host )^ were typically ±0.25 magnitudes in the region 
of the Mb, z plane where the AGN are found, but increased 
to ±0.5 mag where the nucleus was dominant. Similarly, er- 
rors in nuclear magnitudes were ±0.25 mag for bright nuclei 
but as much as ±0.5 mag for host-dominated objects. 

In total, nuclear Mb magnitudes were obtained for 66 
AGN in the sample (10 had no detectable nuclear compo- 
nent), and these data form the basis for the calculation of 
the OLF below. Nuclear absolute magnitudes were found to 
lie in the range —14.6 > Mb > —24.1. The region of the 
AGN Mb{jomc),z plane sampled in this study is shown in 
Fig. 1. 

Since we are attempting to contruct an OLF from an 
X-ray-selected sample we need to ensure that there is a good 
correlation between nuclear optical and X-ray luminosity for 
objects in the SBL sample. We used the monochromatic 
2keV X-ray luminosity, Z/2k e v( n uc) and 2500A UV fluxes, 

* In the SBL study all magnitudes were based on the AB sys- 
tem. Nuclear B(AB) magnitudes were derived by applying a mean 
(B — I) (AB) = 0.2 colour correction to the nuclear /(AB)-band 
magnitudes obtained from the fit to the HST data. For objects of 
this colour, there is a negligible colour term between B and B(AB) 
passbands. For the purposes of this analysis we have therefore as- 
sumed Affl(AB) = Mb for the nuclear regions. 
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Figure 1. Nuclear absolute magnitude v. redshift for the 66 AGN 
with point source detections in SBL. 

^2500A(nuc), listed in the SBL paper. Z/250C>A(nuc) is based 
on the fitted nuclear Ms(nuc), assuming a power-law op- 
tical/UV spectrum of the form /„ oc v~ 5 . Errors on the 
X-ray flux range from 5 per cent for the brightest X-ray 
sources to 25 per cent for the faintest X-ray sources (Gioia 
1990). Although the X-ray luminosity for each source is ex- 
pected to be dominated by the AGN, it is impossible to rule 
out some contribution from the host galaxy, particularly for 
the lowest luminosity sources. 

The least-squares fit to the observed relationship be- 
tween Z/2kev(nuc) and L25ooa(uuc) plotted in Fig. 2 gives 
a relation of the form L2k c v(nuc) oc L25ooA(nuc)°' 82±0 ' 08 , 
consistent with other studies (e.g. Green et al. 1995). 

2.2 The l/V a OLF estimate 

Space densities were derived using the 1/V a estimator 
(Avni & Bahcall 1980). The OLF was constructed from the 
summed contributions of n AGN using: 

n 

$(M B ),z) =J2yi S ( M b ~ Mb) 

i=l a 

Va being the accessible co-moving volume of the i th AGN. 
The estimate of V a was based on z m ax derived from the 
AGN's X-ray flux and EMSS flux limits, assuming an X-ray 
spectral index ax = 1; f v on v~ ax . To construct an optical 
OLF, we binned the 1/V a estimates according to their opti- 
cal nuclear Mb magnitudes. We computed Poisson errors on 

/ vO.S 

the binned estimates of the OLF using a = ( y] J 

Not all 127 AGN with z < 0.15 in the EMSS were ob- 
served by SBL and so a straightforward normalising factor 
of 0.6 (76/127) was applied to the area coverage function 
when computing the accessible volume. KS tests confirmed 
that the redshift and flux distribution for the SBL sample 
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Figure 2. The correlation between nuclear UV and X-ray flux 
from the central component. The line is the slope of best fit calcu- 
lated as 0.82 ±0.08. Typical error on the flux estimates are shown 
in the lower right hand corner. 



Table 1. Nuclear Optical Luminosity Function 



Mb (nucleus) $ (Mpc 3 mag 1 ) Objects 



-24.00 


4.30 ±4.30 x 10~ 8 


1 


-23.00 


3.98 ±3.98 x 10~ 8 


1 


-22.00 


1.31 ± 7.59 x 10~ 8 


3 


-21.00 


5.61 ± 1.91 x 10~ 7 


10 


-20.00 


1.43 ±5.69 x 10~ 7 


14 


-19.00 


1.90 ± 5.63 x 10~ 7 


16 


-18.00 


1.57 ±7.40 x 10~ 7 


7 


-17.00 


3.42 ± 1.63 x 10~ 6 


8 


-16.00 


7.80 ±5.18 x 10~ 7 


3 


-15.00 


3.38 ± 2.54 x 10- 6 


2 


-14.00 


3.36 ±3.36 x 10~ 6 


1 



is consistent with the sample being drawn at random from 
the z < 0.15 EMSS parent sample. 

The OLF was calculated at 1-mag intervals for the full 
redshift range z < 0.15. We made no correction for evolution 
across the redshift bin. We also constructed separate OLFs 
for AGN in elliptical and spiral hosts to investigate any host- 
related trends. 



3 RESULTS 

The differential OLF calculated for an Einstein-de Sitter 
universe in which Ho = 50/iso km s -1 Mpc -1 , Qm = 1, 
J7a = 0, is presented in Table 1 and plotted in Fig. 3, over 
plotted with data from K97 (their Table 5) . As a comparison 
values using total galaxy luminosity (host + nucleus) are 
shown in Table 2 and plotted in Figure 5. 

The K97 data comprises 27 objects extending out to a 



Table 2. Total Optical Luminosity Function 



M B (host ± AGN) $ (Mpc~ 3 mag" 1 ) Objects 



-24.00 


4.30 ±4.30 x 10~ 8 


1 


-23.00 


1.19 ±0.69 x 10~ 8 


3 


-22.00 


9.92 ± 2.55 x 10~ 7 


17 


-21.00 


4.91 ± 1.31 x 10~ 6 


28 


-20.00 


6.98 ±3.69 x 10" 6 


14 


-19.00 


2.48 ±2.30 x 10~ 6 


2 


-18.00 


1.09 ± 1.09 x 10~ 6 


1 



redshift of 0.3. Of these, eight were at redshifts greater than 
the z = 0.15 cut-off adopted in the SBL sample, all of which 
have Mb < —24, i.e. brighter than the most luminous AGN 
in the SBL sample. 

In Fig. 3 we have also plotted two predictions of the 
z < 0.15 OLF based on the luminosity evolution models of 
Boyle et al. (2000) . These authors fit a variety of evolution- 
ary models to a data set comprising over 6000 QSOs with 
M B < -23 and 0.35 < z < 2.3 selected from the 2dF QSO 
redshift survey (Boyle et al. 1999) and the Large Bright 
QSO survey (LBQS, Hewett et al. 1995). Boyle et al. (2000) 
found that luminosity evolution models provided acceptable 
fits to the data, with exponential evolution (L* tx e kT ) as 
a function of look back time (r) favoured for a go = 0.05 
universe and as a general second order polynomial with 
redshift (L* oc in^^*^ 2 ) for q = 0.5. The extrapolated 
z < 0.15 OLFs for the best-fitting 'exponential' and 'poly- 
nomial' models are shown as the short- and the long-dashed 
lines respectively in Fig. 3. The model OLFs have been 
plotted over the magnitude range consistent with the corre- 
sponding range (with respect to M B ) over which they were 
derived at z > 0.35. We obtained a reduced \ 2 = 1-0 for 
the exponential model fit to the SBL data at Mb < —19, 
but were able to reject the extrapolation of the polynomial 
model at the 99 per cent confidence level. 



4 DISCUSSION 

There is good agreement both in slope and normalisation 
between our estimate of the OLF and the K97 OLF at 
Mb < —20. However, at fainter magnitudes the two esti- 
mates diverge. Our estimate of the OLF turns over to a 
much flatter slope whereas the K97 OLF continues to rise 
steeply. However, there are only three AGN in the K97 sam- 
ple with Mb > —20, whereas the SBL sample contains 37 
AGN at these fainter magnitudes. It is therefore most prob- 
able that the difference between the two data-sets (2a) at 
these magnitudes is simply due to small number statistics 
in the K97 sample. 

It is possible that the X-ray selection used to generate 
the SBL dataset is systematically biased against AGN with 
low optical luminosity. However, the correlation between 
L2kev(nuc) and L2500A(nuc) plotted in Fig. 2 demonstrates 
that there is no systematic trend for objects with lower 
optical luminosities to exhibit relatively weaker X-ray-to- 
optical flux ratios. Indeed, the derived relation L 2 k c v(nuc) oc 
L2500A(nuc)°' 82±0 ' 08 implies the reverse, i.e. that AGN with 
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Figure 3. Binned OLF of the 66 AGN in the SBL survey (filled 
circles) compared with data from K97 (triangles). The solid line 
denotes the least squares fit to the data. Predicted z < 0.15 OLFs 
based on luminosity evolution models of Boyle et al. (2000) arc 
also shown (short dashed line: 'exponential' model, long dashed 
line: 'polynomial' evolution, see text for details). 



Figure 4. The relation between nuclear and total galaxy luminos- 
ity for the SBL datasct. At fainter magnitudes the distribution is 
very broad, with no clear correlation between nuclear and total lu- 
minosity. Typical errors for the nuclear component are indicated; 
errors are greater at fainter nuclear magnitudes. 



lower optical luminosities have stronger X-ray-to-optical flux 
ratios. 

In common with other groups (La Franca & Cristiani 
1997, Goldschmidt & Miller 1998), K97 have used their de- 
termination of the low redshift OLF to claim that the slope 
of the bright end of the OLF ($(L)) flattens significantly 
from <E>(L) oc L~ 3 ' 6 at z > 0.6 to $(L) oc L~ 2,5 at z < 0.3. 
Such an observation would rule out pure luminosity evolu- 
tion models, in which the shape of the OLF remains invari- 
ant with redshift. This is in marked constrast with our result 
that an extrapolation of the exponential form of a pure lu- 
minosity evolution model derived at z > 0.35 still provides 
an adequate fit to the z < 0.15 OLF. 

To investigate the discrepancy between these results, 
we fitted our binned estimate of the z < 0.15 OLF with a 
two-power-law model of the form: 

$(L) kT L > V 

oc L L < L* 

Fixing a 'break' luminosity at V = M B — —20.5, we de- 
rived slopes of a = —2.1 ± 0.3 and f3 — —1.1 ± 0.1 using 
a weighted least squares technique. This fit is over-plotted 
as the solid line on Fig. 3. This slope for a is indeed flat- 
ter than that derived at high redshift and is consistent with 
other estimates of the slope of the low redshift OLF, includ- 
ing the most recent determination of the bright end slope of 
the X-ray QSO LF (a = -2.6) by Miyaji et al. (1998). How- 
ever, the value of a derived in this crude fashion is strongly 
dependent the choice of M%, and the inclusion of the two 
brightest bins that each contain a single object. By choosing 
a 'break' magnitude of Mb = —21.5, and restricting consid- 
eration of the data points in the OLF to those bins which 



contain more than one object, we can increase this value to 
a = —2.6 ± 0.3. This is, admittedly, a very crude analysis 
and more sophisticated fitting of a smooth two-power-law 
function similar to that used to fit the higher redshift OLF 
would yield a more accurate estimate of the statistical errors 
associated with fitting the OLF. 

However, it is also likely that systematic errors play 
an equally important role in the determination of the local 
OLF. We attempted to estimate the sizes of such errors by 
first exploring the factor used to correct for the host galaxy 
luminosity. 

In the SBL sample the host galaxy luminosity was ex- 
plicitly removed using fits to the individual HST images. For 
the bulk of their sample (i.e. 0.07 < z < 0.3, or M B < -22) 
K97 relied on a two-step procedure using corrected CCD 
magnitudes measured in an aperture of diameter approxi- 
mately equal to that of the seeing disk, hence replacing total 
magnitudes with small aperture magnitudes. Subsequently, 
a further correction factor for host galaxy luminosity was 
applied by adopting a template host galaxy of Mb = —21. 

In Fig. 4 we have plotted nuclear luminosity against the 
total galaxy luminosity for the SBL sample. As also found by 
della Ceca et al. (1996), although the more powerful AGN 
reside in the more luminous hosts, the distribution of the 
ratio between nuclear and total luminosity is not constant; 
indeed the spread becomes very large at total absolute mag- 
nitudes fainter than ~ Mg(nuc) = —22.5. 

We re-computed our estimate of the z < 0.15 OLF 
based on the SBL dataset using total, instead of nuclear, 
absolute magnitudes. The resulting OLF is shown in Fig. 5, 
with the both the original fit derived for the OLF and the ex- 
trapolated model fits plotted as a comparison. We find that 
although the bright end of the OLF has steepened apprecia- 
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Figure 5. Binned LF of the 66 AGN in the SBL using total 
absolute magnitudes. Solid line: least-squares fit to original OLF; 
Short dashes line: OLF predicted by exponential evolution model; 
long dashes: OLF predicted by polynomial evolution model, see 
text for details. 



bly toa~ —3, both model fits are now clearly incompatible 
with the OLF computed in this fashion. 

The treatment of galaxy luminosities can thus result in 
significant differences to the estimate of the OLF at zero 
redshift. Note also that the evolutionary models have been 
derived from high redshift OLFs uncorrected for host galaxy 
light. Although the assumption that the host galaxy light 
makes an increasingly small contribution to the total lumi- 
nosity of QSOs at high redshift may well be correct, some 
spectacular counter-examples have already been discovered 
(Aretxaga et al. 1995, Brotherton et al. 1999). 

We conclude that the large statisitical and systematic 
errors associated with the determination of the low redshift 
OLF and extrapolation of the OLF at higher redshifts make 
it difficult to rule out luminosity evolution models on the 
basis of shape of the low redshift OLF. 

We also attempted to derive nuclear OLFs for AGN 
with bulge-dominant (E/S0) and disk-dominant (spiral) 
hosts. Following SBL, the distinction between the two broad 
classes and elliptical was made on the basis of the bulge-to- 
total light ratio, B/T, for the host galaxy. Galaxies with 
B/T > 0.5 were classified as E/S0 (44 in the sample), those 
with B/T < 0.5 as spiral (22 in the sample). The nuclear 
OLFs for different types of host galaxy is plotted in Fig. 6 
(after re-normalising the E/SO LF to the same space den- 
sity as the Sa/Sb LF). There are no significant differences 
between the two OLFs; confirmed by a KS test on the cu- 
mulative luminosity distributions for spiral and E/S0 hosts. 

This result is consistent with observation by SBL that 
B/T for the host galaxy was independent of nuclear lumi- 
nosity. In contrast, we note that McLure et al. (1999) found 
that the fraction of elliptical hosts increased significantly 
amongst the highest luminosity AGN. However, these ob- 
servations are still based on relatively small datasets and 
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Figure 6. OLFs for the SBL sample split according to host 
galaxy. The points for the E/SO hosts are displaced to fainter 
magnitudes by 0.1 mag for clarity. 



the McLure et al. (1999) study predominantly samples a 
different luminosity regime (Mb < —23) to that under in- 
vestigation in this analysis. 



5 CONCLUSIONS 

Results from a comprehensive, unbiased X-ray selected sam- 
ple of AGN using the 0.1 arcsec resolving power of the HST, 
have enabled the first direct estimate of the nuclear OLF for 
AGN to be constructed. 

The OLF derived illustrates a two power law form sim- 
ilar to that derived for QSOs at higher redshifts and as 
such is different to the largely featureless power-law OLF 
claimed for the low redshift AGN identified in the Ham- 
burg/ESO QSO survey. However, any discrepancy only oc- 
curs at Mb > —20, where previous estimates of the OLF 
from the Hamburg/ESO survey are dominated by statisti- 
cal errors arising from small number statistics. 

The OLF is consistent with an extrapolation of the ex- 
ponential pure luminosity evolution derived at z > 0.35 
by Boyle et al. (2000), although the 'best-fit' slope for the 
bright-end slope of the OLF is flatter than predicted by such 
pure luminosity evolution models. 

Given the large uncertainties associated with current 
estimates of the low redshift OLF (not least in the present 
analysis) and the extrapolation of evolutionary models to 
low redshift, it is almost certainly premature to rule out lu- 
minosity evolution on the basis of the current determinations 
of the low redshift OLF. 

Further detailed imaging work on optically-selected 
samples of low-moderate redshift AGN/QSOs will clearly 
help resolve this issue. With the superior imaging capabil- 
ity of the new generation of ground-based telescopes (Keck, 
Gemini) we may look forward to such data being obtained 
in the near future. 
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